Modes of magnetic resonance in the spin liquid phase of CS2CUCI4 
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We report the observation of a frequency shift and splitting of the electron spin resonance (ESR) 
mode of the low-dimensional 5* = 1/2 frustrated antiferromagnet CS2CUCI4 in the spin-correlated 
state below the Curie- Weiss temperature 4 K but above the ordering temperature 0.62 K. The shift 
and splitting exhibit strong anisotropy with respect to the direction of the appUed magnetic field 
and do not vanish in zero field. The low-temperature evolution of spin resonance response is a result 
of the specific modification of one-dimensional spinon continuum under the action of the uniform 
Dzyaloshinskii-Moriya interaction (DM) within the spin chains. Parameters of the uniform DM 
interaction are derived from the experiment. 



Quantum antiferromagnets on triangular lattice ex- 
hibit a variety of unusual phases with intriguing physical 
properties. In this work we report results of electron 
spin resonance (ESR) study of quasi two-dimensional 
(2D) triangular antiferromagnet CS2CUCI4 in the para- 
magnetic phase. CS2CUCI4 realizes distorted triangular 
lattice with strong exchange J along the base of the tri- 
angular unit and weaker exchange integral J' — 0.34J 
along lateral sides of the triangle. The material orders 
below Tn — 0.62 K into a long-range 3D spiral state [ij. 
This Ncel point is far below the Curie- Weiss tempera- 
ture Tew K. In the intermediate temperature range 
Tn < T < Tew a spin-liquid phase with strong in-chain 
spin correlations takes place. Magnetic properties of 
CS2CUCI4 are well described by a quasi-one-dimensional 
(ID) model of weakly coupled spin S = 1/2 chains. In 
particular, inelastic neutron scattering experiments [ll, Q 
are nicely explained by a parameter-free theory of dis- 
persing continuum of one-dimensional spinous [3|. More- 
over, the same description, extended in Ref. [j] to in- 
clude residual symmetry-allowed Dzyaloshinskii-Moriya 
(DM) interactions between the Cu-spins, is able to consis- 
tently explain majority of the phases in the T — H phase 
diagram, for three different directions of the magnetic 
field 0. In addition, several numerical investigations 
[1, 0] of the spatially anisotropic J — J' Heisenberg model 
also find that in a wide region of the parameter space 
< J' /J < 0.6 the spin chains are almost decoupled 
and inter-chain spin correlations are exponentially weak. 
This behavior follows from highly frustrated structure of 
the inter-chain exchange J' and provides additional jus- 
tification to the ID-based approach Q to CS2CUCI4. 

Electron spin resonance (ESR) represents a sensitive 
probe of low-energy magnetic excitations with zero mo- 
menta, q « 0. It is well known as a method pro- 
viding fine details of low-energy spin spectra in both 
ordered and paramagnetic [8l4lOl| phases. Importantly, 
both the position and the linewidth of the ESR signal 



are crucially sensitive to anisotropics of the spin Hamil- 
tonian. Here we report a dramatic manifestation of the 
uniform Dzyaloshinskii-Moriya (DM) interaction on the 
ESR spectrum of CS2CUCI4 in the paramagnetic phase. 
We observe and explain splitting of the ESR mode into 
two modes with lowering the temperature, significantly 
extending previous study [llj, and also find strong de- 
pendence of the ESR absorption on the polarization of 
the microwave magnetic field h with respect to the DM 
vector D in zero external field. 

The crystal samples were grown by two methods. First 
set of samples was prepared by the crystallization from 
the melt of a stoichiometric mixture of CsCuCla and CsCl 
in an ampule, which was slowly pulled from the hot (T= 
520°C) area of the oven. The second method was crys- 
tallization from solution [l2l |. The room temperature 
crystal lattice parameters for both sets of crystals are 
a=9.75 A, b=7.607 A, c=12.394A. The samples obtained 
by last method have natural faces and are lengthened 
along b-axes. Both sets of crystals demonstrated identi- 
cal ESR signals. The study in the frequency range 9-90 
GHz was performed with a set of homemade spectrom- 
eters equipped with cryomagnets. The crystal samples 
were placed within microwave resonators at the maxi- 
mum of the microwave magnetic field. The resonance 
lines were taken at the fixed frequency, as dependences of 
microwave power, transmitted through the resonator, on 
the external magnetic field. Temperatures down to 1.3 K 
(0.4 K), were obtained by pumping vapor of "^Hc ("^He), 
correspondingly. Cooling down to 0.1 K was achieved 
with the help of a dilution cryostat Kelvinox-400. 

At temperatures above T = 10 K crystals of CS2CUCI4 
exhibit conventional paramagnetic resonance of Cu^"*" 
S" = 1/2 magnetic ions with a single narrow line cor- 
responding to g-factor values ga = 2.20 ± 0.02, gi, = 
2.08±0.02, gc = 2.30±0.02 for the magnetic field applied 
along a, b and c axes correspondingly. These (7-factor 
values agree well with those reported earlier at T = 300 
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FIG. 1. ESR frequency at T = 1.3 K for H \\ b. Dotted line 
is paramagnetic resonance with g = 2.08, solid line - theory 
(see text). Empty circles stand for resonances loosing inten- 
sity with cooling. Upper insert: gap vs temperature, various 
symbols correspond to different cryostats. Lower insert: evo- 
lution of the resonance line with cooling. 
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FIG. 2. ESR frequencies at T = 1.3 K for H \\ a. Dotted 
line is paramagnetic resonance with g = 2.20, dashed line - 
theory (see text). Upper insert: lineshape at T = 1.3 K, 
V = 21 GHz. Points - experimental data, solid line is the 
result of a fit by the sum of two Lorentzians (dashed lines). 
Lower insert: evolution of the resonance line with cooling. 



K and 77 K [13|, and, for g^, at T = 4.2 K 

Significant evolution of the ESR spectrum was found 
on cooling below 6 K. This evolution depends strongly on 
the microwave frequency v and orientation of the exter- 
nal magnetic field H. (i) For H || b, we observe a single 
Lorentzian line shifting with cooling to lower fields (see 
Fig. [T]). At 1^=14 GHz the resonance field is near zero. 
At lower frequencies, we observe a strong decrease of the 
intensity with cooling, and no visible shift from the para- 
magnetic resonance field, in contrast to ESR lines taken 
at V >17 GHz, which exhibit a shift and do not loose in- 
tensity at cooling. At v—'d GHz the integral intensity at 



T=1.3 K is a half of that at T=4 K. The resonant frequen- 
cies above 14 GHz may be well fitted by the frequency- 
field relation of an ordered antiferromagnet with "gap" 
A/(27r;i) = 14 GHz, 



(1) 



{ii) For H || a, c, the ESR line strongly broadens on 
cooling, and its shape distorts, as shown in the inserts 
of Fig. [2] The absorption at T—1.3 K may be fitted 
by a sum of two Lorentzians, indicating the splitting of 
the resonance mode. The i^(-ff) dependence at H |j a 
is presented in Fig. [21 Two modes at T = 1.3 K were 
resolved for measurements in the frequency range 14-50 
GHz. At higher frequencies the splitting is masked by 
natural linewidth, which prevents the resolution of the 
doublet, while at lower frequencies the signal becomes 
too broad. The angular </> dependence of the resonance 
field Hr = if (sin 0, cos (p, 0) within the a — b plane, taken 
at 1^=27 GHz and T=1.3 K, is shown in Figll For H || c 
the splitting of about 0.35 ±0.1 T is resolved for frequen- 
cies 27 and 31 GHz. 

(Hi) We have also studied polarization dependence of 
the ESR absorption a.t H 0. Fig. 2] shows transmis- 
sion vs field dependences at T — 1.3 K and T = 6 K for 
different orientations of the microwave and static fields 
with respect to the crystal axes, but with h _L H. The 
measurements are done at J/ = 14 GHz. At T > 6 K 
both the position and the lineshape of the paramagnetic 
resonance line are the same for all three principal polar- 
izations of the microwave field. However, at T = 1.3 K, 
the absorption in zero field for h |j b is at least three 
times more intensive than for h |j a, c. 

The splitting of the ESR line for H || c and the po- 
larization dependence of the signal were also detected 
previously in Ref. fu\ . 

We now turn to explanation of the data and show that 
all unusual features of the ESR signal reported above 
are naturally explained by quasi-lD nature of CS2CUCI4. 
The relevant Hamiltonian 

x,y,z 

-g^sH • S:r,y,z + ■•■ (2) 

contains three terms: the first describes intrachain ex- 
change J {x runs along crystal b axis), the second - uni- 
form DM interaction 'Dy^z between chain spins, and the 
third is the usual Zeeman term, allowing for anisotropic 
5-factor. The dots stand for the omitted interchain ex- 
change as well as DM interactions on interchain bonds. 
Detailed symmetry analysis of the allowed DM interac- 
tions Q shows that there are four different orientations 
of the DM vector (see Fig. 6 in Q) depending on chain's 
integer coordinates y,z: r)y.z = Da{—l)^a + Dc{—l)^c. 
Here z indices magnetic b — c layers, while y numerates 
chains within a layer. Crucially, crystal symmetry forbids 
DM vector to have component along the chain b axis.Q 
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FIG. 3. The angular dependence of the resonance field Hi in 
the ab plane for v = 27 GHz at T = 1.3 K. Thick line is the 
theoretical prediction with _Da/(4ft) — 8 and Dc/{4:h) — 11 
GHz. The measured values of the resonance field deep in the 
paramagnetic phase (at T = 10 K) are presented by crosses, 
the dashed line is a theoretical fit of a paramagnetic resonance 

with ga,b,c of CS2CUCI4. 
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FIG. 4. The polarization effect at i/ = 14.44 GHz for H ^ 0. 
Circles with error bars show the uncertainty resulting from 
normalizing the signals to zero absorption a.t H > 1.5 T. 



The essence of the observed ESR line splitting can be 
explained by considering a single Heisenberg chain with 
uniform DM interaction with vector D along the local 
z axis, when the field is lined up along the same axis 



Dz, H — >■ Hz). In this geometry the DM inter- 



action can be gauged away by position-dependent rota- 
tion of lattice spins. 



The rotation angle tan(a) 



0+ iax qz 
x,y,z ' '^x.y.: 



(3) 



x,y,z- 

-D/J is chosen so 



as to eliminate the DM coupling from the Hamiltonian. 
The transformed % describes spin chain with quadratic 
in D easy-plane anisotropy. To linear in D/J accuracy 
the ESR response of this model coincides with that of 
an isotropic Heisenberg chain, detailed study of which is 
described in Ref . [loj . ESR absorption is determined by 
the transverse structure factor 5a;a;(a;,q) = Syy{uj,q) of 
the chain, evaluated at g = 0. This is produced by char- 
acteristic to the spin-1/2 chain spinon continuum and is 
described by a sum of delta-function peaks at frequen- 
cies huj± = \hvq±g^j:BH\, see eq.(3.16) of [l^l and FiglS] 
Here v — TTjao/{2h) is the zero-field spinon velocity and 
flo is the lattice spacing |14| . The response in the original 
(un-rotated) basis is obtained by un-doing the momen- 
tum boost (O, which corresponds to setting q = D/Jao 
in the expression for the structure factor S^x/yy This 
immediately implies the splitting of the ESR line into 
two lines, at frequencies hoj± = {guBH ± TrD/2\. Note 
that for small D/H the splitting ttD = tt{DI + Dlf/'^ is 
linear in D which justifies our neglect of the easy-plane 
anisotropy. 

While the four-sublattice structure of the DM vec- 
tors makes it impossible to study the D || H configu- 
ration in CS2CUCI4, the above argument allows one to 
immediately understand polarization-dependent absorp- 
tion (finding {Hi) above) in zero magnetic field, at the fre- 
quency 27ri/ = TrD/2h. Since the microwave absorption is 
proportional to the square of the microwave field compo- 
nent perpendicular to the effective field, we conclude that 
configuration with h _L D (that is, h || b) should result in 
maximal possible absorption. For h along the a (c) axis, 
absorption is a factor D^/D^ {D^/D^) smaller. Using 
numerical estimates of Da,c derived below, the absorp- 
tion for different polarizations should obey the relations: 
Pb/Pa — 1-6 and PbjPc — 2.8. This agrees qualitatively 
with our data presented in Fig. 2]- the zero field absorp- 
tion at = 14 and 17 GHz is a factor of 3 more intensive 
for h II b than for h || a, c . 

To understand our findings (i) and (ii), one needs to 
analyze the arbitrary orientation of the H and D vectors, 
which is the problem solved in [lij . The main physi- 
cal point of Ref.[l5| consists in the observation that a 
uniform DM interaction, much like a spin-orbit Rashba 
interaction in 2D conductors [ll-[ll, acts on spinon exci- 
tations of ID chains as an internal momentum-dependent 
magnetic field. Technically, this observation follows from 
spin-current formulation [lj| of the problem ([2]) , which is 
valid below the strong-coupling temperature scale To ^ 
Je-"'"^ {S = 1/2 for the spin chain here)^. The 
spin currents M^/^ represent spin density fiuctuations 
of spinons near the right/left (R/L) Fermi points of a 
ID system. It follows that the right and left-moving 
currents experience different, in magnitude and direc- 
tion, total magnetic fields: B/j = H -I- hvD / gjisJ and 
Bl = H — hvT) / gjisJ ■ It is then natural that the ESR 
response of such a chain consists of two peaks at frcquen- 
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FIG. 5. Spinon spectrum of S" = 1/2 Heisenberg chain for q ~ 
with (sohd hnes) and without (dashed lines) DM interaction. 
When D II H, momentum is boosted by g = D/ Jao. 



cics 27rfti'R/L = 5Ms-Bfl/L.[15l 

For the Cs2CuCl4-specific configuration of chain DM 
interactions and H = (Ha, Hb, He) our analysis predicts 

{2nhiyn)^ = [gblJiBHbf + {gaPsHa + (-l)"^i?a/2)' + 
{gc^iBH, + {-^)y^^D,/2f, (4) 

{2TTh,^l^f = {gtiiBHbf + {gafisHa - {-iyTTDa/2f + 
(<7,MBffc-(-l)%i?c/2)^ (5) 

Observe that while vr ^ vl for fixed chain indices {y, z), 
these chiral ESR frequencies transform into each other 
under lattice translations )• (y±l,z±l). These 

equations naturally explain the difference between H || 
a, c and H || b situations. For H || b the external and 
the internal DM field are mutually perpendicular and the 
vector sum H ± D has the the same absolute value for 
both signs. One then finds single ESR frequency = 
of the form ([1]) where the gap is in fact determined by 
DM interaction, A = ^T^/Dl + D^/2. 

For H II a, c the DM field has a component along H 
and the frequencies of the two spin excitations are differ- 
ent. In principle, because of the four values of the Dj, ^ 
vector, four different frequencies may be present when the 
field H is oriented within the a — c plane. We attempted 
to check this experimentally but failed to resolve sepa- 
rate lines within the broad band of absorption. 1J| For 
H within the a — b plane ^ and ([S]) predict two dif- 
ferent frequencies. Solving equations (|4I5P at the fixed 
frequency for the magnetic field we were able to fit the 
experimental angular dependence in Fig. [3l In this way 
we obtained Da/ (Ah) = 8±2 and £'c/(4/i) = 11±2 GHz, 
which corresponds to 0.29 ± 0.07 and 0.39 ± 0.07 T. We 
have to note here that for H || c the observed splitting of 
the frequency doublet Ai/ = 0.35 T is about 50% smaller 
than the prediction of ()4I5|) with these Da/c values. This 
deviation from the predicted value is dicussed in 14 1 . 

A more detailed comparison between the experiment 
and the theory would require extending the latter to fi- 
nite temperatures since the reported measurements are 
performed in the intermediate temperature range Tn = 
0.62K < T < Tew = 4K. Temperature dependence of 
the gap A for H |j b in a wider range including T/v is 
shown on FiglTJ The variation of the gap in the range 



lK<T<2Kis much slower than right near Tn. This 
supports our main assumption that the reported spectral 
features of CS2CUCI4 reflect specific spin chain physics 
which dominates the temperature range Tn < T < Tew- 

In conclusion, we demonstrated that the uniform DM 
interaction, which is a distinctive feature of CS2CUCI4, re- 
sults in a new kind of spin resonance in a S= 1/2 chain an- 
tiferromagnet. The observed spectrum is a consequence 
of the splitting of the chain's spinon continuum by the 
internal magnetic field which is produced by the uni- 
form DM interaction. We believe that a similar phe- 
nomenon is possible in a higher-dimensional magnetic 
systems with fractionalized spinon excitations. Our find- 
ings differ strongly from the well-known single-frequency 
resonance of = 1/2 antiferromagnetic spin chain with 
staggered DM interaction [l^l . It also differs from a con- 
ventional 3D antiferromagnet which acquires an energy 
gap in the spectrum only below the ordering transition. 
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Supporting material for 
"Modes of magnetic resonance in the spin liquid 
phase of CS2CUCI4." 

In this work we report results of electron spin reso- 
nance (ESR) study of quasi two-dimensional (2D) an- 
tiferromagnet CS2CUCI4 in the paramagnetic phase. 
CS2CUCI4 realizes distorted triangular lattice structure, 
shown in Figll] with the strong exchange J along the 
base of the triangular cell and with the weaker exchange 
integral J' for the exchange paths along lateral sides of 
the triangle, their ratio is J'/J = 0.34. In-chain DM 
interaction 

T>y^,^Da{-ira + D,{-l)yc, 

derived in Ref. [l|, describes four possible orientations of 
the D vector on the chain with integer indices y, z, see 
Figure [1] Here y and z axes are aligned along a and 
c crystallographic axis of an orthorombic Pnma lattice, 
respectively. 
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FIG. 1. Parameters of spin interactions in CS2CUCI4. Four 
different orientations of tire chain DM vector D are illustrated 
by red arrows. 



ANGULAR DEPENDENCES OF THE 
RESONANCE FIELDS 

In addition to the angular dependence of the resonance 
field on the angle </) within the ab-plane, presented in the 
main text, we present in Fig 15] the angular dependence of 
the resonance field on the angle within the the bc- 
plane. In this case the agreement is only qualitative: 
the splitting at H || c, obtained by a two-Lorentzian 
fit constitutes about 50% of the predicted by the DM- 
theory. Nonetheless, the predicted splitting of the reso- 
nance fields is within the broad band of absorption, ob- 
served at H II c, as indicated by blue arrows in the upper 
panel of Fig. [2l 

Fig IS] presents the evolution of the absorption band for 
the field oriented within the ac-plane. According to our 
theory (see main text) the ESR line should split into 
four resonances, see FigH] We failed to resolve four lines 
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FIG. 2. Upper panel: v = 21 GHz resonance lines for the 
magnetic field H in the bc-plane. The curves are taken at 
the equidistant angles. Blue arrows for H || c curve indi- 
cate the calculated resonance fields for the splitted ESR line. 
Lower panel: The tp angular dependence of the resonance 
field in bc-plane for v ^ 27 GHz at T = 1.3 K. Thick line 
is the theoretical calculation (see text) with Da /(4ft) = 8 
and Dc/{4:h) — 11 GHz. Dashed line shows calculated angu- 
lar dependence for a paramagnet with g- factors of CS2CUCI4, 
crosses represent experimental values at T = 10 K (deep in 
the paramagnetic phase). 

because the width of the individual lines is too broad. 
Nevertheless, the expected positions of four lines, marked 
by arrows in Fig|3] do lie within the measured absorption 
bands. 



FREQUENCY-FIELD AND TEMPERATURE 
DEPENDENCES FOR H | c 

The frequency- field dependence for H || c at T = 1.3K 
is shown in the Fig IS] At T = 1.3 K the splitting is 
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resolved only for two frequencies in the middle of the 
frequency range. At low frequencies, the large linewidth 
masks the splitting, while at higher frequencies the split- 
ting predicted by the DM-theory is not observed at this 
temperature. However, at lower temperatures the split- 
ting of the ESR line for H || c becomes more visible. It is 
clearly seen even at the high frequency edge of the range, 
v = 76 GHz, as is seen in Fig. [6l The splitting observed 
at T = 0.79 K equals 0.6 ± 0.1 T, which corresponds sat- 
isfactory to the predicted value of 0.75 T (horizontal shift 
between the asymptotic parts of dashed lines in FiglS]). 
For H II a and at T = 0.8 K, we observe similar splitting 
of 0.6 ± 0.1 T for the frequency u = 72.9 GHz. 



TEMPERATURE DEPENDENCE OF ESR LINES 
AT FREQUENCIES ABOVE AND BELOW THE 
GAP A FOR H II b 

It was mentioned in the main text that at the frequen- 
cies below the DM gap, i.e. v < A = M GHz, the ESR 
line docs not exhibit a shift but instead looses the inten- 
sity. This observation is illustrated in Fig. [3 where the 
data for v = 9.63 GHz ESR signal are presented. One 
can clearly see that the resonance field does not change 
with lowering the temperature, while the intensity of the 
line does diminish strongly. 

To contrast this with quite different evolution of the 
signal at higher frequencies > A, we present in Figl7]de- 
tailed temperature evolution of the ESR line aX v = 26.9 
GHz. Here the shift of the line's center to the lower val- 
ues of the resonance magnetic field is obvious. The signal 
also broadens such that the overall intensity is conserved. 

For other field orientations, in particular for H || a, c, 
the ESR line completely disappears for low frequencies 
V < A. 
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FIG. 3. The ESR lines for v = 27 GHz and T = 1.3 K at 
different orientations of the static field within the ac-plane. 
Arrows mark the resonances, expected at T = for decoupled 
spin chains according to the DM theory (see main text) with 
Da/ Ah = 8 and _Dc/4fi = 11 GHz. 



THEORETICAL DETAILS 
Lattice Hamiltonian for D || H 

When H || D, unitary rotation (3) of the paper trans- 
forms chain Hamiltonian (2) into 

X 

~gtJ^BH~si, (1) 

where for brevity we do not write chain's integer co- 
ordinates y, z. As claimed in the paper the obtained 
Hamiltonian ([1]) possesses an easy-plane anisotropy of the 
strength D^/2 J. Neglecting it takes us to the Hciscnberg 
limit discussed in the paper and illustrated in Fig. 5 there. 



Low-energy theory 

The obtained Hamiltonian ([T|) can also be used to de- 
rive the low-energy form of the DM coupling. This is 
done by bosonizing ^ (with H = 0) which results in Q 

H^Jdx ^{{d^^)^ + {dj)^}, (2) 

where ((f, 9) is the pair of conjugated bosonic fields which 
describes low-energy excitations in the rotated basis. In 
particular, the transverse component of the spin density 
operator (at momentum tt) is expressed via ^ as 0| 

S+ = A^e'P^ (3) 

where /3 = 27ri? determines scaling dimension of the spin 
operator in terms of the known compactification radius 
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FIG. 4. Resonance fields for p = 27 GHz ESR as a function 
of the angle in the ac-plane according to the DM-theory with 
Da/4:h — 8 and Dc/4:h = 11 GHz. Dashed line shows calcu- 
lated angular dependence for a paramagnet with g-factors of 
CS2CUGI4, crosses represent experimental values at T = 10 K 
(deep in the paramagnetic phase) 



R, see Refi, and A3 is non-universal amplitude. Also 
note that spinon velocity w in ([2]) is function of magnetic 
field. It monotonically decreases with H and reaches zero 
at the saturation 0]. Near H = the v{H) dependence 
is weak and we can safely approximate v by its zero-field 
value everywhere in this study. Thus, when necessary, 
we use 



2~ 



(4) 



for the velocity v. Here is the lattice constant in units 
of which we measure coordinate x. 

Comparing ^ with equation (3) of the paper we im- 
mediately observe that at low energies that unitary ro- 
tation simply corresponds to the shift of momentum by 
a = —D/J (which is exactly the physics discussed in the 
main paper). Hence bosonic field 9 of the original theory 
and 9 of the rotated one are simply related as 



D 



= u —X. 



/3J 



(5) 



Hence the low-energy theory of the original Heisenberg 
chain with uniform DM interaction can be obtained from 
(HD via (P so that 



n 



dx -m^f + {dj 



D 

Jj 



I 



dx -m^f 



?} 

vD 



{djy} + -T-fdxO + const. (6) 



f3J 



Using that at if = /3 = V27r and M| 



-d^9/p 



we find that the last term in (jB]) is proportional to the 



difference of the right and left spin currents. This allows 
us to identify the lattice and the low-energy forms of the 
DM interaction as 



• X S 



x+l 



vD 



dx (A/^ - Ml) (7) 



This result is written in terms of fluctuating spin densi- 
ties (spin currents [1]) M'^^^ defined in terms of (Dirac) 
fermions ^r/l near the right/left (R/L) Fermi points, 
correspondingly. More precisely, 



B./L;s'- 



(8) 



It turns out that the spin-current formulation is partic- 
ularly convenient for the problem at hand (see [Hj for 
pedagogical explanation of the technique) . 

To complete the low-energy description we also need 
the low-energy form of the Zeeman coupling, which is 
well-known and reads (for the moment we keep magnetic 
field along z axis as well) 



V2^ 



dx dx^p 



= g^iBH / dx {M'j^ + Ml). 



(9) 



In writing ([2]) we have omitted a number of marginal 
terms which account for the easy-plane anisotropy as well 
as residual backscattering interaction between right- and 
left- spin-currents 0, @1 • These terms do not affect our 
main result - the splitting of the ESR signal into two 
due to the interplay of the external H and the internal 
DM D fields. They are, most likely, important for under- 
standing the width of the individual ESR line, which is 
the problem of much greater theoretical complexity and 
is left for future studies. 

With the help of ([7]) and ^ we can now write down 
the low-energy Hamiltonian of the Heisenberg chain sub- 
ject to both Zeeman and DM interactions acting along 
different directions. 



chain 



Ho - y dx gfiBH{M^ 



M'l) + 



dx {vD/J){M'^ 



Mt) 



(10) 



where Hamiltonian "Ho of the ideal Heisenberg chain is 
written in Sugawara form Q as 



-Ho = — / dx {Mr ■ Mr + Ml -Ml} 



(11) 



To connect with previous expressions, it is worth noting 
that is in fact equivalent to the abelian bosoniza- 
tion form ([2]). The more complicated- looking (jlip has 
one significant advantage over it makes evident 

a remarkable emerging SU(2)i^x SU(2)l symmetry H, 
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which allows for independent rotations of the right- and 
left- spin currents. This key for the following symmetry 
emerges below the strong-coupling energy /temperature 
scale To Je'''^ [S = 1/2 for the spin chain here)0, 
when the spin current formulation becomes appropri- 
ate. Above that temperature a more conventional semi- 
classical description of the spin- 1/2 chain in terms of 
fluctuating Neel vector order parameter field is valid . 
In fact, one can think of Tq scale as of 'Haldane' tem- 
perature below which the difference between integer and 
half-integer spin chains becomes pronounced. It is the 
temperature (or, cquivalently, energy) scale below which 
the true spinon nature of the quantum ground state of 
the critical spin- 1/2 chain is manifest. 

Note that upper indices d and h on currents M/j/^ 
in (fTUl) denote directions of the vectors D and H. Ob- 
serve that DM coupling (last term in "Hchain) is odd under 
spatial inversion which interchanges R and L - this sym- 
metry distinction is already evident in the lattice Hamil- 
tonian (2) of the main paper, see also 



90 
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FIG. 5. ESR frequencies at T = 1.3 K for H \\ c. Dotted line 
is paramagnetic resonance with g = 2.30, dashed line - DM 
theory (see text), Da/4fi = 8 and D^f-ih = 11 GHz. 



Hchain =Ho- giiB / dx{BRKPj, + BiAf|) 



ESR as a chiral probe 

The remaining steps closely follows Ref. 0. We now 
perform independent rotations of the right- and left- cur- 
rents in (fTOl). so that in the rotated basis 



(12) 



the Hamiltonian describes unusual situation whereby 
right and left-moving currents experience different^ both 
in magnitude and direction, magnetic fields: Bj^ — 
H + hvD/gfiBJ and = H — hvD/gfiBJ- Observe 
that such a rotation does not affect Ho which retains its 
quadratic in spin currents form (jlip . 

It is now rather natural (see Refs. Q and H for de- 
tails) that the ESR response the spin chain consists of 
two peaks at frequencies (here we make use of ^) 



^irh fR/L = 5MB|Bfl/i| Lg^sH ± 7rD/2|. 



(13) 



The obtained result includes H || D and H = situations 
discusses in the main paper as special cases. 

It is interesting to note that in the described case of a 
single spin chain with constant D vector the ESR exper- 
iment becomes chiral probe: (|13|) shows that right- and 
left- moving excitations can be accessed independently by 
simply varying (and rotating) external magnetic field H 
with respect to DM axis D. In the case of CS2CUCI4 this 
interesting point is lost as the DM axis takes on four dif- 
ferent directions, depending on the y, z coordinates of the 
spin chain. As a result the right ESR frequency, equation 
(4) of the main paper, for, say, chain with even coordi- 
nates y, z coincides with the left frequency (equation (5)) 
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FIG. V — 75.88 GHz ESR lines at different temperatures 
for H II c. 



from chains with odd y, z indices. Thus, the crystal struc- 
ture of the material masks the chiral nature of the ESR 
probe. 



Frequency dependence of the ESR response 

We would like to describe here one possible reason for 
a pronounced frequency dependence of the ESR signal 
as described in Section „ above. Consider first H || b 
arrangement when our analysis predicts that the signal 
is centered about the single frequency 



2T:h VR,L - ^/{^l^^HYT^^^niTDf)^. (14) 
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FIG. 7. Upper panel; set of 9.63 GHz ESR lines at diflferent 
temperatures for H \\ b. The narrow resonance at ~ 0.35 T 
is the signal of g=2.0 DPPH marker; a tiny resonance, ap- 
pearing only at low temperatures, is due to a paramagnetic 
impurities. Lower panel: set of 26.92 GHz ESR lines at differ- 
ent temperatures for H \\ b. The narrow resonance at ~ 0.9 
T is the DPPH signal. 

This is just equation (1) of the main paper (see also dis- 
cussion following equations (4,5) there). This equation 
clearly states that there is a minimal ESR frequency, 
given by i^b.min = A/(27r?i) = D/{4h), below which 
no resonance is possible in this geometry. Experimen- 
tally the resonance is observed both above and below 
i4),min, but with markedly different temperature depen- 
dencies. The low-frequency < I'b.min) response re- 
mains paramagnetic-like and quickly looses intensity on 
lowering T, see Figure [T] On the contrary, the high- 
frequency response (Figure [7]) follows prediction (fT4|) 
rather nicely, as detailed in the main text. 

We'd like to argue that the observed difference is a 
finite-T effect. Eq. ([TUl describes idealized T = sit- 
uation of completely decoupled spin chains, to which 
CS2CUCI4 is only an approximation in the intermedi- 
ate temperature window between Tq discussed above and 
the ordering temperature T/v . In this temperature inter- 
val one expects sizable thermal population of spinous at 



pretty much all momenta, not restricted to the q = 
limit in which the current field theory is formulated. 
These thermal spinons arc not particularly correlated 
and can be thought of as forming an ideal gas of neu- 
tral fermions. It seems physically sensible to think that 
their ESR response is paramagnetic and not sensitive to 
the low-energy effects due to the internal DM field. It is 
this response that is measured sX v < fb,mm in Fig. [T] 
According to this argument the intensity of the param- 
agnetic part of the signal should diminish with lowering 
T, simply following diminishing of the population of ther- 
mal spinons, in agreement with our observation in Fig.[7l 
The higher-frequency response {v > t'b.min) is instead the 
intrinsic property of the spin chain. 

It is easy to extend this argument to other field direc- 
tions. Let, for example, H || a: then according to equa- 
tion (5) of the paper we can achieve the situation when 
gafJ-sHa = TTDa/2 and 2TThvp>{z = odd) = 2TThvL{z = 
even) = -kDc/^^ while 2'Khh'fj{z = even) = 2'KtivL[z = 
odd) = tt-^/D^ + 4D2/2. Thus in this case z^a.min = 
Dc/ (4/i) and this condition is realized for 1 /2 of chains 
in the system. Similar argument shows that for H j| c 
the minimal resonance frequency is t'c,min = £'a/(4ft). 
To reiterate, below the corresponding minimal resonance 
frequency we expect standard paramagnetic ESR signal 
due to thermally excited population of spinons, and the 
intensity of such a signal should decrease with lowering 
of the temperature. 

It is interesting to note that there is one 'magic' ori- 
entation of the external field for which fmin = 0. This 
happens for gaPsHa = ■nDa/2,H}, = and gcHsHc = 
ttDc/S, in which case 1/4 of the chains experiences zero 
total field. It would be interesting to study this special 
orientation in more details in the future. 
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